NASA/TM— 20 10-216368 



Conservational PDF Equations of Turbulence 


Tsan-Hsing Shih 

Ohio Aerospace Institute, Brook Park, Ohio 
Nan-Suey Liu 

Glenn Research Center, Cleveland, Ohio 


July 2010 


NASA STI Program ... in Profile 


Since its founding, NASA has been dedicated to the 
advancement of aeronautics and space science. The 
NASA Scientific and Technical Information (STI) 
program plays a key part in helping NASA maintain 
this important role. 

The NASA STI Program operates under the auspices 
of the Agency Chief Information Officer. It collects, 
organizes, provides for archiving, and disseminates 
NASA’s STI. The NASA STI program provides access 
to the NASA Aeronautics and Space Database and 
its public interface, the NASA Technical Reports 
Server, thus providing one of the largest collections 
of aeronautical and space science STI in the world. 
Results are published in both non-NASA channels 
and by NASA in the NASA STI Report Series, which 
includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant phase 
of research that present the results of NASA 
programs and include extensive data or theoretical 
analysis. Includes compilations of significant 
scientific and technical data and information 
deemed to be of continuing reference value. 
NASA counterpart of peer-reviewed formal 
professional papers but has less stringent 
limitations on manuscript length and extent of 
graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies that 
contain minimal annotation. Does not contain 
extensive analysis. 

• CONTRACTOR REPORT Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, often 
concerned with subjects having substantial 
public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific and 
technical material pertinent to NASA’s mission. 

Specialized services also include creating custom 

thesauri, building customized databases, organizing 

and publishing research results. 

For more information about the NASA STI 

program, see the following: 

• Access the NASA STI program home page at 
http://www.sti.nasa.gov 

• E-mail your question via the Internet to help@ 
sti.nasa.gov 

• Fax your question to the NASA STI Help Desk 
at 443-757-5803 

• Telephone the NASA STI Help Desk at 
443-757-5802 

• Write to: 

NASA Center for AeroSpace Information (CASI) 
7115 Standard Drive 
Hanover, MD 21076-1320 


NASA/TM— 20 10-216368 



Conservational PDF Equations of Turbulence 


Tsan-Hsing Shih 

Ohio Aerospace Institute, Brook Park, Ohio 
Nan-Suey Liu 

Glenn Research Center, Cleveland, Ohio 


National Aeronautics and 
Space Administration 


Glenn Research Center 
Cleveland, Ohio 44135 


July 2010 


Acknowledgments 


This work was supported by the Subsonic Fixed Wing Project as part of the NASA Fundamental Aeronautics Program. 
The authors highly appreciate the fruitful discussions with Dr. P. Tsao. 


This report is a formal draft or working 
paper, intended to solicit comments and 
ideas from a technical peer group. 


This report contains preliminary findings, 
subject to revision as analysis proceeds. 


Trade names and trademarks are used in this report for identification 
only. Their usage does not constitute an official endorsement, 
either expressed or implied, by the National Aeronautics and 
Space Administration. 


This work was sponsored by the Fundamental Aeronautics Program 
at the NASA Glenn Research Center. 


Level of Review. This material has been technically reviewed by technical management. 


Available from 


NASA Center for Aerospace Information 
7115 Standard Drive 
Hanover, MD 21076-1320 


National Technical Information Service 
5301 Shawnee Road 
Alexandria, VA22312 


Available electronically at http://gltrs.grc.nasa.gov 


Contents 


Abstract 1 

1.0 Introduction 1 

2.0 PDF of Turbulent Variables 2 

2. 1 PDF of Turbulent Velocity fu(V; x,t) 2 

2.2 PDF of Turbulent Species x,t) 3 

2.3 Other Turbulent Quantities and Functions of Turbulent Variables 4 

2.3.1 Statistical Mean of Analytical Function of Turbulent Variable 5(0,) 4 

2.3.2 Statistical Mean of Derivatives of Turbulent Variable dUt/dxj and d 2 Ui/dxjdxj 4 

2.4 Joint PDF 5 

2.4. 1 Joint PDF and Marginal PDF 5 

2.4.2 Conditional PDF and Conditional Mean 6 

2.4.3 Relationship Between Unconditional Mean and Conditional Mean 7 

2.4.4 Summary 9 

3.0 Transport Equation for Turbulent Velocity PDF V; x,t) 10 

3. 1 Mean and Differentiation Commutation 1 0 

3.2 Taking Mean on Navier-Stokes Equations 1 1 

3.3 Traditional Velocity PDF Equation 12 

3.4 Conservational Velocity PDF Equation 13 

4.0 Transport Equation for Turbulent Species PDF x,t) 15 

4.1 Taking Mean on Species Equation 15 

4.2 Traditional Species PDF Equation 17 

4.3 Conservational Joint Species-Velocity PDF/j^PA)/; x,t) Equation 17 

5.0 Concluding Remarks 19 

Appendix A. — Fine Grained PDF, Velocity and Species PDF Equations 21 

A.l Transport Equation for Fine Grained PDF f'(V; x,t) 21 

A.2 Transport Equation for Velocity PDF fu{V; x,t) 21 

A.3 Transport Equation for Joint Species-Velocity PDF/, ,®(E,V|/; x,t) 23 

Appendix B. — Miscellaneous Formulations 25 

References 26 


iii 


NASA/TM— 20 10-216368 




Conservational PDF Equations of Turbulence 


Tsan-Hsing Shih 
Ohio Aerospace Institute 
Brook Park, Ohio 44142 

Nan-Suey Liu 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


Abstract 

Recently we have revisited the traditional probability density function (PDF) equations for the 
velocity and species in turbulent incompressible flows. They are all unclosed due to the appearance of 
various conditional means which are modeled empirically. However, we have observed that it is possible 
to establish a closed velocity PDF equation and a closed joint velocity and species PDF equation through 
conditions derived from the integral form of the Navier-Stokes equations. Although, in theory, the 
resulted PDF equations are neither general nor unique, they nevertheless lead to the exact transport 
equations for the first moment as well as all higher order moments. We refer these PDF equations as the 
conservational PDF equations. This observation is worth further exploration for its validity and CFD 
application. 

1.0 Introduction 

To explore the turbulence modeling employed in the PDF method, we revisited the constructions of 
the traditional PDF equations for turbulent velocity and species, which were described by several 
researchers in great details, for example, by Pope (Ref. 1). In those constructions, a Delta function (which 
is referred as the fine grained PDF) was extensively used, together with a few impressive mathematical 
techniques to invoke the Navier-Stoke equations into the identity relationship developed from the Delta 
function. Since the entry point for the Navier-Stokes equations was through the “conditional mean”, two 
or three conditional means appeared in the velocity PDF equation and additional two conditional means 
appeared in the species PDF equation. In these traditional PDF equations, all the conditional means were 
considered as unknowns and modeled empirically. 

In the present study, we started out by taking the statistical mean directly on the Navier-Stokes 
equations to form an integral equation, from which a PDF equation is then constructed. Depending on the 
way of taking mean on some particular terms, especially the pressure gradient VP and the molecular 
diffusions vV 2 Uj , TV 2 0 ; , the resulted PDF equation can end up with different forms. Briefly speaking, 
if we treat, for example, the molecular diffusion term v V 2 (7, as a separate random variable from U„ then 
its mean (vV 2 Uj^ will be expressed as an integration of the product of the velocity PDF and the 
conditional mean, i.e., 


(vV 2 !/,-) = J V (v 2 f/ / |F)/(L ; x,t)dV , 

where the conditional mean (v V 2 L'/| V 'j is a new u nk nown. However, if we take the mean based on the 
definition and the commutation rule between mean and differentiation, then we have 


’V 2 U i ) = vV 2 (U i ) = vV 2 \vj{V ;xj)dV = ;x,t)dV , 
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which is in closed form. The former way will lead to the traditional PDF equations (see Sections 3.3 and 
4.2). The latter way will lead to the closed equations for the velocity PDF (see Section 3.4) and the joint 
species-velocity PDF (see Section 4.3). Since a sufficient but not necessary condition was involved in 
their construction, the resulted PDF equations are neither general nor unique. Nevertheless, we observed 
that they can lead to the exact transport equations for not only the first moments (i.e. mean velocity and 
species, etc.) but also all higher order moments deduced from the Navier-Stokes equations. This feature is 
referred as the conservation of the PDF equations. 

Later, we also found that, for the traditional way of utilizing the “fine grained PDF”, all the 
conditional means (e.g.,/v V 2 U t \ , ^vV 2 0,-|\|/^ , etc.) can also be systematically modeled through 
similar sufficient condition, hence the traditional PDF equations can be closed exactly as the PDF 
equations proposed in the present study. The details are described in the Appendix A. 

Section 2.0 introduces the basic definition of probability density function of turbulent velocity, 
species, other turbulent quantities and their joint PDF; the relationships between the marginal PDF, the 
joint PDF and the conditional PDF; and the relationship between the mean and the conditional mean. 
Section 3.0 demonstrates how to construct the velocity PDF equation starting directly from the Navier- 
Stokes equations. Section 4.0 demonstrates the construction of the transport equations for the species PDF 
and the joint species-velocity PDF. 

2.0 PDF of Turbulent Variables 

2.1 PDF of Turbulent Velocity fu(V\ x,t) 

Turbulent velocity L) (x ; - , t) is a random vector variable, its probability density function PDF at a 
single point (in space and time) is denoted as fjj (V; x,t) = fu (V\,V 2 ,Vt, ; x\,X 2 ,xt, , t ) , where 
V =(V\,V 2 ,V^) is the sample space coordinates of U„ and x = (xi,x 2 ,* 3 ), 1 are the location coordinates 
and the time. Formally, fjjiV', x,t) is a scalar function of V and x,t, but its arguments have different 
physical meaning and mathematical transformation. The V represents a point in the velocity sample space 
( - x < Vj < go, i = 1, 2, 3 ) at which the probability density f { < of U, is defined. The x,t indicate that f v is a 
field and a process. It is important to note that V is independent of x and t. The probability density 
function fu (V ; x, t ) itself is not random and is fully defined at each V for every single point x,t, it is a 
differentiable (hence continuous) function of V, x, t. Furthermore, its integration over the whole sample 
space V at every single point x,t must be equal to one, because the total probability for all events must be 
100 percent, i.e., 


GO GO GO 

J J j f u (V 1 ,V 2 ,V 3 ;x 1 ,X 2 ,x 3 ,t)dV 1 dV 2 dV 3 =l. (1) 

—GO —OO —OO 


Or written as 

\f v (V;x,t)dV = l. 


( 2 ) 


Where J ( 


)dV is an abbreviation of 


J 




) dVidV 2 dVi • It is also important to note that the 


—CO — OO — OO 

argument V in f L , represents ‘ V \ , V 2 , V 3 ’ so that fu(V; x,t) is a j oint probability density function 
.A/ ( Ij , L 2 , k) 1 -V, I ) of the turbulent velocity components U\,U 2 ,U$ . 
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As Pope (Ref. 1 ) pointed out that, with the PDF fjj (V ; .v, t ) , we can define various one-point 
statistical properties of the random velocity U, (x, I ) , for example, the mean velocity (or the first moment) 



(U i ) = \v i f u (V;x,t)dV 


( 3 ) 


All other higher order moments at any single point of x,t (also referred as one point moments or one 
point correlations) can also be defined: 


(U'Vj^jViVjfuiV; x,t)dV 

(UjU jU k ) = \v t Vj V k fu ( V ; x,t)dv (4) 


In particular, the Reynolds stresses liijUj ') defined by the Reynolds decomposition u { = U t ~(Ui) are 
fully defined by /[/(F; x,t) : 

{•‘.Uj)=((u l -{u l ))(uj-(uj)))=(u l uj-u l (u 1 )-u i {u l ) + (u l ){u j )) =(u l Uj)-(U l ){Uj) (5) 

Equation (5) indicates that the turbulent kinetic energy k - (u i Uj')/2 is fully defined by the velocity PDF 
fu(V; x,t). 


2.2 PDF of Turbulent Species x,t) 

Turbulent scalar variables, for example, the species O, (x, t), i = 1,2,-- • n , are random variables. The 
PDF of species is denoted as /$ (v| /; x,t) , where v| / = \\ij is the n-dimensional sample space coordinates of 
the speciesO,. Therefore, /®(\j/; x,t) is a joint scalar PDF of multi species. Analogous to the turbulent 
velocity, the one -point statistical properties of the turbulent species ®, are fully defined by the PDF 
/®(y]/; x,t) . For example, the mean is written as 

(®, ) = J\|/ ; -/®(\)/;xT)</\|/, i' = 1, 2, •••,/? (6) 

And the higher order moment as: 

(®i) = JV/ /®(V; *,t)dy 

(® 1 ® 2 ...® n ) = | V |/ l¥2 --.v|/„/ (1) (\|/; x,t)d\\i 


In particular, the i th species fluctuating intensity ? 'j defined by the Reynolds decomposition 
4* / = ®, - (® ; ) is fully defined by /®(v|/; x,t) , 

)(*>!) ( 8 ) 
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2.3 Other Turbulent Quantities and Functions of Turbulent Variables 

There are several other turbulent quantities, for example, the pressure P, the chemical reaction source 
S(®,), the pressure gradient dP/dXj , the velocity gradient VC/,- , and the higher order derivatives such as 
V 2 Uj , V 2 ®, . They may be viewed as other separate random (scalar or vector) variables, in addition to 
the turbulent velocity C/, and the species ®,. We will see later that these “random variables” often appear 
as a “conditional mean” in the traditional PDF equations. 

There are two types of quantities: one is the analytical function of the turbulent variables such as 
S(®,) = ®“®? • • • , and the other is the derivatives of the turbulent variables such as V 2 ®,- . The quantities 
S(®,) and V 2 ®, are both random because ®, is random. From the probability theory the statistical 
properties of S(®,) and ®, can be determined by the same PDF /®(v|/i X P) ■ Flowever, the statistical 
properties of V 2 ® ; is not that straightforward. 

2.3.1 Statistical Mean of Analytical Function of Turbulent Variable A(®,) 

The statistical mean of S(®,) is defined as 

(SX®,)) = J A(v|fi)/ct>(v|/; x,t)d\/ (9) 


where S(\) /,) is the functional form in the sample space V|/ for S(®,). The function S( ) is the assumed 
analytical function. 

2.3.2 Statistical Mean of Derivatives of Turbulent Variable dUjldxj and ('fU l ldjxpx J 

There are at least two ways to take means on the derivatives of velocity. Let us first use the 
commutation rule (Section 3.1) to express the following means: 

( VC/, ) = V (U, ) = V J V, f v (F; () iV - J V, 
x ' J dxfdxf 


Lienee, the mean of derivatives can be determined by the derivatives of the PDF of turbulent variable. 
Equation (10) can be further expressed as 


<™.> 

(V2[/,) = J 


Vi dfu 

fu dx 


fu(V; x,t)dV 


J J 


Vi d 2 fg 
fu dxjdxj j 


fu(V]x,t)dV 


or 


{VU i ) = \Dj{V i )f u {V-x,t)dV 
(v 2 U i ) = \L(V i )f u (V;x,t)dV 


where D/ V t ) and L( V,) are defined as 
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( 12 ) 


D m= 

Wd = 


Vi dfu(V ; x,t) 
fu(V;x,t ) &,■ 

Vi d 2 fu(V;x,t ) 


fu ( V ; jc, f) dxjdxj 


Analogy to Equation (9), we may view D,( F,) and L( K,) as the functional form s of VC/,- and V 2 C/,- in the 
sample space V. 

Another way to take mean on the derivatives of velocity is to view them as some different random 
variables from the velocity; hence it will involve the joint PDF or the condition means of the velocity 
derivatives. This will be described in the next Section, where the relationships between the conditional 
means and the functions Df V,j, L(V , ,) will be revealed. 


2.4 Joint PDF 

2.4.1 Joint PDF and Marginal PDF 

The Navier-Stokes equations contain several random variables, for example, the velocity the 
pressure P (or the pressure gradient VP) and the molecular diffusion vV 2 f, . The scalar transport 
equations contain the velocity U,, the species ®, and its molecular diffusion TV 2 ®,- . Hence a joint PDF 
for the joint variables of U„ ®„ P, V 2 Uj =L U ‘ , V 2 ®,- = , i.e. f U Q> P pj L <s> (F,\| t,p, l u ,/® ; x,t), 

may be needed in the analysis together with the marginal PDF, such as fjj (V ; x, t) , (\|/; x, t ) , and the 

joint PDF of fu <d (F,V| /; x,t) . The arguments p, l u , are the sample space variables of P, V 2 Uj and 
V 2 ®,- , respectively. 

By definition, PDF must satisfy the normalization condition: 

\\\\\fu,<b,P,L u ,ip(y ’V’P’l 17 ; x,t)dVd\\idp dl u dl® =1 (13) 

The relationship between the joint PDF fuc&ppJ (F,\| r,p, l u ,/° ; x,t) and the marginal PDF 

fu(V;x,t),f 0 (\| /; x,t), f P (p; x, t ) , f L u (l u ; jt, t ) f L ® (/° ; jt, t ) are defined as 

fu(V;x,t) = ffflA p L u IP (V,\\f,p, l u , /®; x,t)d\\idpdl u dl <s> 

/o(v; x ,t) = IWa o P,L U , l? iV,^f,p, l u ,/® ; x,t)dVdpdl u dl® 

f P (p;x,t) -JfJR a> p L u l? (Vi'tyiP) l U J®; x,t)dVd\ydl u dl® (14) 

fju (l u ;x,t) = | J J J fu,<s>,p,i¥ (V, \|/, p, l u , /® ; x,t)dVdydpdl ® 
fpt, (l®;x,t) = J J J J fu,<b,p,L u ,IP (V’V'P’ l U ,l®;x,t)dVdx\idpdl u 

This is clearly required by the normalization condition for any type of PDF, including the marginal and 
joint PDF. With this basic concept, we can obtain various relationships between different joint PDF, for 
example, the joint PDF of U-, and ®, can be defined as 

fuM V ’V’ x ’ t ) = ,<S,P,L U ,iP (F,\| !,p, l u ,/°; x,t)dpdl u dl ® (15) 

It should be noted from Equations (13) and (14) that the statistical mean of a random variable can be 
equally define by either the joint PDF or its marginal PDF, for example, 
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( 16 ) 


(Ui ) = 1 1| 1 1 Vjfu ^ pju L <b (V,\\r,p,l u ,1®; x, t)dVd\\rdp dl u dl® 

= \Vi[ y \\\\fu,<s>,P,LU J* ( v > V| j,p,l u ,l 9 ;x, t)d\\rdp dl u dl® ) dV 
~ jVj fj(V ; x,t)dV 

(®i) = J J J J J x \’ifu,<bJ\lP r ,L« l u ,1®', x,t)dVdy/dpdl u dl® 

-MWK <t> P IP I? l u ,1® x,t) dVdpdl u dl®^d\\f (17) 

= jV//®(v;*>0<fy 

'•* mJh V jfu O P IP Z* (^V'76 l u \ x, t)dVd\\idp dl u dl® 

= \\ Vi V* (III fu,o,P,L u jP (V’V’P’ X,t)dpdl u dl^dVdx y (18) 

= W v i^ifu,<b(y,\f\x,t)dVd^ 

Equations (16), (17) and (18) indicate that taking mean of a random variable only needs the relevant 
marginal PDF or the relevant joint PDF. 

2.4.2 Conditional PDF and Conditional Mean 

Fet us examine an expression: 

\Vifu^{V,Mr,x,t)dV (19) 

which is taking mean on the random variable (7, at a fixed Vj/, i.e. O, = V|/ ; - . If we define the ratio of the 
joint PDF to the marginal PDF as 


fu\<t><y v; x,t) 


fuM V ’V’ x ^ 

fb(W> x d) 


( 20 ) 


Then the expression (19) can be written as 

| v i fu,o (^ Vi 0 dV = (\|/; x, t) • | Vj f v ^ (V\ v|/; x, t) dV t 

= /<d(v; x,t)-(Ui\i/f) 


( 21 ) 


Where (Uj \j/^> is an abbreviation of (U,- 1 0 = \j/^> and is referred to as the conditional mean of U, on the 
condition of ®, = \| /,• : 


{Ui^lVifu^ty^xtfdV (22) 

And fjj |(j,(F|v|/; x,t) is referred to as the conditional PDF of (7, on the condition of ®, = vjq , which 
satisfies the normalization condition because of the definition (20): 

\fu\^{V\ T ,x,t)dV = \ (23) 
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Similarly, the conditional mean of the random pressure P(x,t ) on the condition of U, = V} is denoted as 
(P(x,t)\v) or (/ 7, | l 7 ) and the corresponding conditional PDF is denoted as 
/® P L u z® U (V- P' \V , ; x,t ) , then we may naturally write 


(P(x,t) \V) = JJJJ'W. L <s,p(\\r,p, l u ,/® \V , ; x,t)d\y dpdl u dl® 

= \pfp\u(p\v ; x,t)dp 


(24) 


The above conditional PDF is defined as 




fu,^>p,ip,ip w ’V’P’ > x d) 

fu(V ; x,t) 


(25) 


or 


f PU (p v ;x,t)= 


fuj><y>P\x,t) 
fu(V ;x,t) 


(26) 


In Equations (24) and (25), the sample space variables V|/, l u , /® are irrelevant. 

2,4.3 Relationship Between Unconditional Mean and Conditional Mean 

It should be noted from Equation (24) that a conditional mean is a function of the sample space 
variable that has been excluded from the process of conditional mean (e.g., Vj in this case). Obviously, 
(P\V) is a function of sample space variable Vj , and taking ‘mean’ on (P\v) over the sample space Vj 
will lead to its unconditional mean (or just mean) pressure ( P ) , which is a function of Xj, t .: 

\{ F \ v )fu<y\x>*)dV = KW'-i u(p\V ; x,t)dp^fu(y ; x,t)dV t 

= J J P fu,p (V,p;x,t) dVdp 
= ( p ) 


Or vice versa, 


(p) = \\pfuA V ’P ; x,t)dVdp = \fu(y ; x,t){P{x,t)\v)dV (27) 

Equation (27) is a general formula for two or more different random variables (P and (7, in this case) 
involved in the conditional mean of one random variable (P in this case). 

Now, following Equation (27), we may straightforwardly express the means of the derivatives of 
turbulent velocity VU,- , V 2 U, as 


(vu,) = J/[/(K; x, t)(VUj\v)dV 

(28 

(' V 2 Uj) = \fu(V ; x,t)(v 2 Uj\v)dV 

Comparing Equation (28) with Equation (1 1), a sufficient (but not necessary) condition for this identity is 
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Equation (29) indicates the relationship between the conditional means and functions Dj(Vj ) , L(V,) . It is 
noted here that these equations can also be derived by using the fine grained PDF through similar 
sufficient condition (see Appendix A: Equations (100) and (103)). 

For the mean of pressure gradient ( i dP(x i ,t)/dx j ) , following Equation (27), we may write 



dP(Xj,t ) 
dxj 


V )dV 


(30) 


In the past, the conditional mean of the pressure gradient in Equation (30) has been considered as an 
u nk nown quantity even if P = P (U) was assumed. This conditional mean, like the other conditional mean 
(v 2 u t \v) , is actually carried by the traditional velocity PDF equation and is modeled empirically. Now, 
if we use the rule of commutation and take the mean on the pressure gradient (assuming that 
P(x,t) = P(U(x,t)), or P( x, t ) = P(U( x, t), x, t ) ), we may write 



£-(p) = -^-jp(V)f u (V;x,t)dV 


dx, 


\^-{P(V)f u (V;x,t))d v 


r 1 dP(V)fu 
\fu dxi J 
\P?<y)MV\ x,t) dV 


fu(V;x,t)dV 


(31) 


Here, P(V) is dependent on x,t in general and PP ( V ) is defined as 


Jf(V) = 


1 d[P(Y)fu] 

fu dxi 


(32) 


Comparing Equations (30) and (31), and using a similar argument that leads to Equation (29), we obtain a 
relationship between the conditional mean and the function Ff ( V ) : 


/ dP(Xj,t) 
\ dx t 



' fu S Xi 


(33) 


It is noted here that this equation can also be derived by using the fine grained PDF through the similar 
sufficient condition (see Appendix A: Eq. (104)) 

Now let us look at the mean of the product [/,• -0„ , i.e. (Uj 0 ;? ) . By definition we may write 
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{Ui ® w ) = J ¥« ({ Vifu,<i>(V,y; x,t)dV )t/v|/ 
= J Vn f<b (v; X, t ) (Ui | V|/} c/v|/ 


( 34 ) 


or 


(^; <*>» ) = { v i ( j ¥«/[/,< D (F,\|/; x,t)d\^dV 
~ \ Vj fu(V> x,t)(® n \y) dV 


(35) 


Equations (34) and (35) illustrate that, when taking joint mean, the joint PDF must be used to avoid the 
appearance of the conditional mean. The conditional mean ^(7 ; |v|/^ is also carried by the current 
traditional species PDF equation. 

Finally, let us examine the mean and the conditional mean of a function of the variables O,, say 
iS(®j0 2 • • - ®«) j on the condition of ® = V|/ . Let us first view S(®i®2 ‘ • - ®«) as a separate random 
variable with the sample space variable s which is independent of \|/, then, following Equation (27), we 
write the mean as 


(5(®!® 2 • ■ •<&»)) = \u (v; *> 0 — <&„ )| V W (36) 

The conditional mean in (36) is 

= X,t)ds 

= ^(Vi V2 •••¥«) J/j|® (s\y, X, t) ds (37) 

= 'S'(VlV2 •••¥«) 

Then Equation (36) will end up an expression that is just the definition of its mean: 

(5 , (® 1 ®2 ' ‘ ) = j S(\iiV 2 ■ ■ -Vii)/©(V; x,t)d\f (38) 


This term is related to the chemical reaction rate in the traditional species PDF equation, which is in a 
closed form. Also note that the reason for the conditional mean of S (® 1 ®2 ‘ ' - ® w ) in Equation (36) 
becoming free from the constraint is that 1 S , (® 1 ® 2 ' "® ;i ) is a known function of ®, without involving 
spatial differentiation. 

2.4.4 Summary 

2.4.4. 1 Formulations With the Conditional Mean 


y 2 U i ) = \f u {V-x,t) (v 2 u\v)dv 

(39) 

(VP) = j/c/(F; x,t) (VP\v)dV 

(40) 

v2 ®/) = J/o(¥ ; x,t) (v 2 ® ; \|/)r/\)/ 

(41) 
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( 43 ) 


(5 , i (0i0 2 ---0„) = J/ o (\(/;x,0 (<S)|i|/)<ii|/ (42) 

{Ui ®n ) = J Vr, /© (v; *, 0 (t// 1 v) 

The conditional means (except the one in Eq. (42)) are all considered as unknowns in the traditional PDF 
equations. 

2.4.4.2 Formulations Without the Conditional Mean 

(y 2 u i)=\Vi v 2 fu (v ; x,t) dv 

or (44) 

(v 2 Ui)= \L(V i )f u (V;x,t)dV, L(Vi) = Vi ^ fu 

or (45) 

( v2 0/) = J\)/f V 2 / 0 (v|/; x,t)d\\i 

or (46) 

(v 2 ®,-) = |z,(\|/i)y®(v; x,t) d\\i, z,(\|/ f ) = J 

(,S)(®i®2---®„) = J U (y; x, t) S i (\\i l \\i 2 ---\\i n )d\\i (47) 

(Uf O n ) = J J Vi v) info# (V, v|/; *, t)dVd\f (48) 

These equations are exact and have no u nk nowns other than the PDFs fu, ./<j> and fu$> . The functions 
L(Vj) , L(\\t /) and ff (V) arc all well defined by the corresponding PDF, and the function P( V) is 
assumed as known and will be defined later (Section 3.4, Equation (68)). 


3.0 Transport Equation for Turbulent Velocity PDF f L { V: x,t) 

3.1 Mean and Differentiation Commutation 

It is important to note the commutation property of taking mean and differentiation, that is 

M 20 /A 020 /_2L\ = 2iO 

\dt / dt \dxj j dxj ’ \dxjdxi / dxjdxj’ 


(49) 
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These commutation properties have been shown in several text books, for example by Pope (Ref. 1 ), 
Tennekes and Lumley (Ref. 2). For example, 


dU(x,t) \_/ Um U(x + Axi,t)-U(x,t) \_ Um (U(x + Axj,t))-(U(x,t)) _ d(U(x,t)) 

dxj / >o A Xj / Axj — >0 Axj dxj 


( 50 ) 


where U = {U\, Ui, C/ 3 ) . 


3.2 Taking Mean on Navier-Stokes Equations 

Apply the mean operation on the Navier-Stokes equations for incompressible flow, i.e. 

lw + »ML m _ i^ +vv2(/ \. 

y dt dxj p dxj j 


(51) 


Using the commutation rule, Equation (49), each term in Equation (51) can be expressed as follows, 



d(Ui) _d 


’ 1 = -[v i fu(V;x,t)dV 
dt ot J 




~fu^V ; x,t) 
dt 


dV 


Note, Vj is independent of x- t and t 


(52) 


Similarly, 





v j~fu(V;x,t) 


dV 


(53) 


Taking mean on the pressure gradient term can proceed in two different ways. One is with the conditional 
mean, i.e., via Equation (40): 



Note that the following term in Equation (54) is zero (Pope (Ref. 1) or see Appendix B): 
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The other way without the appearance of the conditional mean is via Equation (45) 


]_dP 

p dx , 


])’-' p \ p ? (y)MV '- x ' ,) W] dv 

= -'j/ r (K 'Ji- [p;(¥y v (V.x.n)dV (56) 




dV : 


Pf(V) 


fu(V; x,t) 


dV 


Note that (see Appendix B): \——(Vi Pf{V) fjj(V\x,t)]dV - 

p J SVj v J ' 


= 0. 


Finally, taking mean on the molecular diffusion term, one way with the conditional mean is via 
Equation (39), 


vV 2 f/ ; -) = J f v (V ;x,t) {vV 2 U t \v)dV = \fu<y',x, t) (v W 2 U 


dV,- 
,\v ) — L 
71 ' dV, 


dV 


= |— [VifuiV ■xd)(vV 2 U j \v))dV -\Vi— [fu<y\ x,t)lytV*Uj\v))dV (57) 


=-\ v ‘ -§r{fv<y ■.*.»{ w ^| f )) 


dV 


Note that (see Appendix B): j" -^-{Vifu(V\ Ar,t)^vV 2 [/y|F^c/F = 0. 
The other way without the conditional mean is via Equation (44), 

dV, ■ 


y V 2 Uij = v j L(Fj)fu(V ; x,t) 


dVj 


-dV 


vl^ViLiV^fu^x^dV-vjv—lLiV^fu^x^dV (58) 


d 


= -j. 


dVj 


(vL{Vj)fu{V-x, tj) 


dV 


Note that (see Appendix B): v f — — (ViL(Vi) x,t))dV = 0. 

J dVj y ’ 


3.3 Traditional Velocity PDF Equation 

Now substituting Equations (52), (53), (54), (57) into Equation (51) and collecting the terms factored 
by Vj , a sufficient but not necessary condition for satisfying the integral equation is that the factored, total 
integrand is zero. Thus, a transport equation for the turbulent velocity PDF fu(V; x,t) is constructed: 
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(59) 


dfu + y dfu _ 8 


dt 


dx; 9Vj 


", / 1 8P 

A) 


V 

t \ P dxf 

1) 


~(fu{^ 2 U,\ v )) 


If the pressure term is decomposed into the mean and fluctuating parts, P - ( P''j + P 

10«P) + P')| 


\_8P_ 

P dxi 


V ) = [ 


dxj 


v ) = -Pal + l 

p dx t \ p dxj 


(60) 


Then Equation (59) can be written as 


VU | ySfu_ lS{P)df v d 
dt 1 dxt p dxi dVi dVi 


fu 


’V 2 Ui V)- 


1 8P' 

p dxi 


(61) 


Equation (61) is exactly the same as the traditional velocity PDF equation (Pope (Ref. 1)). The mean 
pressure in Equation (6 1 ) is considered as a known quantity through its Poisson equation. However, the 
last two terms with the conditional mean were considered as u nk nown and modeled empirically. Using 
the following identity relations (Pope (Ref. 1)): 


, JdP\ 
fu(V; x,t)(— | 




+ ■ 


5 


8V, 


MV; x,t){P' 


8U ; 


8x: 


(62) 


fu(V; x,t)(v 2 U i\v 


dv j L 


Equation (61) can be rewritten as 


a 2 


J dV.dVk 


fu(V ; x,t) 


8Uj 8U k 


8X; 8Xj 


(63) 


«r + ^Mr =vV 2 /t , + i£T>^ + ^i_ 

dt dx, p dxj 8Vj dxidVj 


" Ip' 

Vi 

fu ( — 

V 

l \p 

1 


+ ■ 


1 8 2 


2 dVidVj 


fu 


P_ 

P 


r dUj 8U ^ 


+ ' 

dx-. dx 


J J 


V )- fu 2v 


8Ui dUj 


8x k 8x k 


(64) 


Now we have three unknown terms in Equation (64) that were modeled empirically (Pope (Ref. 1)). 


3.4 Conservational Velocity PDF Equation 

Alternatively, collecting the terms factored by V, from (52), (53), (56) and (58), and following the 
same procedure that leads to Equation (61), we obtain a closed velocity PDF transport equation: 


d fu + y dfu _ 
dt 1 dxj 


8V, 




8V: 


Pf(V ) 

Ju 

P 


(65) 


Here L(Vj) and Pf (V) are known functions: 
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( 66 ) 


Wi) = 


V J y2 fu , p* in 1 d \ p ^fu] 

J fu dxj 


fu JU — / 

In order to determine P(V ) , we first take mean on the Poisson equation 


i V 2 

p dxjdxj 


or 


l^ (f) - gW 

p dxfixj 


By the definition of the mean, we may write Equation (67) as 


- V 2 J P(V)fu (V; x, 0 dV = - f Vj fu (V; x, t ) dV 

L) (JJCj (JJC j 




(67) 


or 


if V 2 (P(V)fu (V; x,t))dV = -\Vi Vj d2 fu(ri*,t) dv 
p J J J dxfiXj 

By equating the integrands on both sides (again, this is a sufficient condition for satisfying the above 
integral equation), we obtain a model equation for P(V) : 


-V 2 [P(nfu] = -V i Vj-f^- (68) 

P CXf Ox j 


Therefore we have constructed a velocity PDF Equation (65) that does not need extra turbulence 
modeling. 

Using Equation (65), it is easy to show that the exact equations for the first moment and higher order 
moments of the Navier-Stokes equations can be deduced by multiplying the Equation (65) with V, and 
Vt Vj respectively and integrating over the velocity sample space V. For example, 


\ v iVj\— + V k ^L 
J | dt dx k 


5 Vi 


{vL(V k )fv)- 


P k S (V) 


dV k 


fu 


>dV 


The left two terms become 

d(U,Uj) | d(u/jfj k ) 

dt dx k 


and the first term on the right hand side is, 

. a . a . Qy. y . 

-\ v i v j—OW k )f u )dV = -\—(v i VjVL(V k )f u )dV + \vL(V k )f u ^-PdV 
OVk O 'k 

= \(vL(Vj)Vifu +vL(Vf)V jfu ) dV 
= v(u i V 2 Uj+U j V 2 U i ) 
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Note that (see Appendix B): — (Vj V.v L(V k ) fu)dV = 0. 

J dV k v ’ 

Similarly, the second term on the right hand side is 


Vi V, 


1 J dV, 


pRv) 


fu 


j 


l TT dP JT 8P ' 

dV — — ( U i f U / — 

p \ dxj 1 Ox, I 


Therefore we obtain the exact equation for the second moment (UjU j j : 

dp' 


d(UjUj) d(u i u j u k ) _ i/ tt dP 

dt dx k 


= -p 


4.0 Transport Equation for Turbulent Species PDF /®(\|/; x,t) 

4.1 Taking Mean on Species Equation 

Taking mean on the transport equation for species O, 


SO,- dd>iU ; . 

+ TV 2 ® / +5 / (® 1 ® 2 ---®„ 


dt dx 


The first term on the L.H.S. of (70) can be written as 


) 


-f 

dtf 


Jv/ 

5 . 

— /©(v; *>0 

dt 


d\\i 

Equation (71) involves a marginal PDF /$ . We can also write this term using the joint PDF fu 

~^T ~ = J* J dvdxv 

= jJVi — (/[/,(D (V,v;x,t)) dVdy 

The second term on the L.H.S. of (70) can also be formed in two ways: 

d(®iUi) d 


dxj 


dx 


Vi 


- J J 'ViVjfuM v >V> x,t)dVdx\i 


d 

dXj 


■/©(v; x,t)(u / K|/ 


d\\r 


or 
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d(®iUj) d 


d x ; 


dx 


lb 


J J ViVjfu , o (V, V; A:, 0 c/Kc/v)/ 


*>0 


( 74 ) 


JF(7v|/ 


The next term V 2 ® ( of (70) can again be written in two different forms: Using Equation (41), we have 

r ( v2 ®i)=rJ/©(v; *,0^v 2 ® ; -|y^<7i|/ = rJ/ 3 )(\|/; ;M)(v 2 ®y|vi^|^-</v|/ 

C\\l j 

_ o . 

= r Wi>(v;*,0(v 2 ® y |v|/) (/(!,( \|/;a:,o(v 2 ® 7 ' | (75) 


cty/ 


= -jV* ^-(/®(v; x,o( r v 2 o y |v)) 


c/V|/ 


Q - - 

Note that (see Appendix B): [ T/./a> (v; x > 0 ( V 2 ® ,• \|/ ) — 0 . 

J dy j L ' 7 1 / J 

Or we can straightforwardly write its mean as 

r ( v 2 ®, ) = r v 2 (® , } = r v 2 { {v,/^ (K, x, o dvdx v 
= j J V/ [ r v2 /(/,<t) (V, V|/; x, o] dVdxsf 
= ||n(v i)fu, o(Y,\\r, x,t)dVd\\i 

Where T(\|/, ) is defined as 

V/ v2 /(/,d> _ V; v2 /d> 


UVi) ■ 


Equation (76) can be further manipulated as 


fu,<b 


/<D 


T (V 2 ® ; ) = T J J L(y j)fu,<& (Y,W> x,t)-^-dVd\\r 

= r J J ^-(v» H^j)fu,o(y,T> x,t))dVd\\f 
-r J JV/ ^-( L(Vj)fu,<b(y>W> x,t)j)dVd\ \i 


=- r Jf 


Vi 


01]/ 


-(L(Vj)fu,<b(y,V> x ’0) 


dVdy 


(76) 


(77) 


(78) 


Note that (see Appendix B): r ll^ -{ViL(\\rj)f Ut ®(V,y; x ,t))dVd\\f = 


= 0. 
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The reacting source term is readily available via Equation (42), 


(Sj (0t0 2 ■ • • )) = J/© ( v; X, 0 (Sj I V|/) d\\i = | /© ( y; x, 0 | \\i) 

= { (79) 


di|/. 


= -J 


¥« 


ch|/ 


-(/o ( V|/; O'S'y (¥iV2 ■ ■ •¥» )) 


(7v)/ 


Note that (see Appendix B): J — — -^ 4 ///(d(v|/; x,0^<S/|v|/^<7y = 0. 
Another expression for the reaction source term can be obtained as 

<5,- (O iO 2 • ■ ■ O „ ) ) = 1 1 5,- (\|/i\|/ • • • \|/„ )/f/ d, (F, v|/; x,t)dVdy 


n 


Vi 


OV(/ 


- ( Sj (\|/f y • • ■ v|/„ )fu,<t> (V, V|/; x,tj) 


dVd\\i 


(80) 


Note that (sea Appendix B):JJ ^(v^tviV-VOA,.^. V. x.»)dVdv - 0 


4.2 Traditional Species PDF Equation 

Collecting terms factored by \|/,- from Equations (71), (73), (75) and (79), and following the 
procedure similar to the construction of the velocity PDF equation, we obtain the species PDF equation: 

f L +^(/»(u>)) = -A-(/ <1 ,((rv2®,|v> + s y (w---v,,))) (si) 

By using Reynolds decomposition Uj = U t - (Uj ) , Equation (8 1 ) becomes 

^ + ^ >l$7 <"'l v » = “ « r v2a, dM') + ■ V,,))) (82) 

Equation (82) is also exactly the same as the traditional species PDF equation derived by Pope (Ref. 1 ), 
which contains two conditional means from the convection term and the molecular diffusion term. These 
terms were considered as unknowns and were modeled empirically. 


4.3 Conservational Joint Species- Velocity PDF _/(;,©( K,\|/; x,t) Equation 

Alternatively, by collecting terms factored by \|/ ( - from Equations (72), (74), (78), and (80) , we obtain 
the following equation: 


dfu ,< D | y dfu , o 
dt j dxj 




5v | )j v ' 


(83) 
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Equation (83) is in closed form. They can be reduced to the marginal species PDF /® (\|/; x,t) equation 
by integrating over the velocity sample space V t . The result is 


MsL+iTj\<h- 

dt ' dxj 



4 (rt( ^ )/ ® ) 


d 

dyj 


(fa>Sj(ViV—Vnj) 


(84) 


In Equation (84), there is an unavoidable conditional mean due to the convection term that involves both 
the velocity and the species. The conditional mean of fluctuating velocity conditioned on O, = \|q- , 

(uj | \|/ ^ , is unknown and must be modeled. 

Now let us return to the subject of joint PDF. Equation (83) is derived from the species diffusion 
Equation (70), the two source terms on the right hand side represent the molecular mixing of “species” 
and the chemistry reaction, respectively. We also have to consider the contribution from the Navier- 
Stokes equations, which can be obtained in a way similar to that leads to Equation (65): 


Pfu,® + y dfu, d> 
dt 7 dxj 


dVj 


(v L(Vj)fu ,®)- 


Pf(V) 


dVj | 


fu, <D 


(85) 


The terms on the right hand side of Equation (85) represent the molecular mixing of “momentum” and the 
pressure source. Therefore, the joint species-velocity PDF fy.^. ( V L , vp ; x, , t ) equation should include all 
source terms on the right hand side of Equations (83) and (85), i.e. 


Qfujb | y dfu, <P 
dt 7 dxj 




a 

8Vj 

d 

dyj 


(pf(V) ] 

JU,® 




( 86 ) 


It is easy to verify that the joint PDF fu <& Equation (86) can be reduced to the equations for the marginal 
velocity PDF fy and the marginal species PDF /$ by integrating over the corresponding sample space, 
and the results are the same as the Equations (65) and (84). It is also easy to verify that Equation (86) can 
deduce the exact equations for the first moments (JJ / ) , (O,- ) and higher order moments (0,0 ^, (U^j 'j , 
etc. 

The functions involved in Equation (86), i.e., L(Vj), P'j’ (V ) , L(\\ij) are all known and given by 
Equations (66), (68) and (77): 


L(Vi) = 


Vj v 2 fu ,® Vj v^fu 


Pf(V) = 


fu,® fu 

1 8[P(V)fu,®~] _ 1 d[P(y)fu] 

fu 


fu,® 


dXj 


dxj 


Wj Pfu,* 

P V 2 fu,® dXjdXj 


ViVj d 2 fu 


(87) 


V 2 /t/ dxidxj 


AV/) = 


yF 2 fu,® _ ViV 2 f® 


fu,® 


f® 
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The Following two relations, which are similar to Equation (12), are also useful: 


d jW) 

Dj(Vi) = 


Vi dfu# _ Vi dfu 


f .<!» OXj f dxj 

Vi d fu,® _ Vi g/o 
fu,<b dx j /<d dxj 


( 88 ) 


The solution of Equation (86) should be able to calculate all one-point statistics for the velocity and 
species, for example, (C/,), (O,-), (UiUj^j, (t/*®,-), (0j0 2 • • •} , etc. We can also calculate other turbulent 

/dO„ 50,. 


, dUi dUj 

quantities, such as ' 


dx k dx k / \ dx k dx k 


(PUi ) and (PSy j , which are related to the dissipation rate 


Ey , s, S(D , the pressure transport and the pressure strain correlations, etc. For example, 


£ i/= v 




'dU ; dUi 


'duA/du,- 


I dx k dx k / \ dx k / \ dx k 


; djc k dx k 


'50„ \ /SO. 


5x/. 


5X; 


k I J 


and 


< W ) - (PUi )~(P) (Ut ) , {psij ) = (PSy ) - (P) (Sy ) 

I Dk(V ‘ )Dk(V J )fu ^ dVdy=\D k {Vi)D k {Vj)fu dV 

= J A- (¥« ) A (¥« )/<d 

(F A ) = J J AA AAc/,0 = J P(V) Vjfu dV 

PSy ) = | J J (P(V)(Dj(Vi) + Di(Vj))f Ut<s ,)dVd V = \\ ( AA( Dy^) + D^Vj^f^dV 


(89) 


(90) 


5.0 Concluding Remarks 

We have constructed a set of conservational PDF equations directly from the Navier-Stokes equations 
and the species diffusion equations through the use of some sufficient but not necessary conditions. 
Therefore, in theory, they are neither general nor unique. Flowever, all these PDF equations are in closed 
form. They are able to deduce the exact transport equations for the first moment and all higher order 
moments. This feature has not been observed from any other existing modeled PDF equations. For 
example, the modeled traditional PDF equations can deduce the first moment equation correctly, but not 
for the higher order moments. 

In this study, we have also defined a few functions, for example, Dj(Vj), L(V,) , Dj (®,) , /.(d), ) 
and Pf (V) , they may be viewed as the models for the conditional means of V U t , V 2 Uj , Vd>, , V 2 ®,- 
and VP . These functions may be used to calculate other turbulent quantities of interest. 
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Future work includes further testing of the present conservational PDF equations for their application 
in the area of CFD and the development of conservational fdtered density function (FDF) equations for 
compressible turbulent flows. 
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Appendix A. — Fine Grained PDF, Velocity and Species PDF Equations 

A.l Transport Equation for Fine Grained PDF f\V; x,t) 

According to Pope (Ref. 1), the fine grained PDF f'(V; x,t) is defined as 


f\V-x,t) = b(u{x,t)-V) 

Taking the following operations, 

df df 8U ) df 8U, 


dt dUj dt dVj dt 

8ff = 8ff8u L ^_8ffdu^^ 

dxj dUj dxj 8 Vj dxj 

We may form an identity equation for the fine grained PDF: 


0 ( , dU t A 
1 dt 


dv iK 

8 r 

~ 8 V, 


f 


dxj j 


d lL +U 

dt J 8x 1 


8f_ 

x j 


8 _ 

8 V, 


f 8U, TT 8U ' 

y dt J 8xj J 


Inserting the Navier-Stokes equation into Equation (93), we obtain a transport equation for the fine 
grained PDF : 


(91) 


(92) 


(93) 


d f\ + u V 


dt 


J dxj " TT 1 


of 


1 8P d 2 Uj 

+ v — 

p dxj 8xjdxj 


(94) 


A.2 Transport Equation for Velocity PDF fu(V; x,t) 

Now, taking mean on Equation (94) we will obtain a transport equation for the velocity PDF. First 
taking mean on the left hand side of (94), 

(er tu <r\ e J /) , d if' u j) 

( dt j dxj j dt dxj 

-V)fu{V'\x,t)dV\^]w' -V)V)fu{V'-x,t)dV' (95) 

_ 8fu(V;x,t) | y 8fu(V;x,t) 
dt j dxj 

To explain the way of taking mean on the right hand side of (94), let us first take the mean on 
fdUi/dXj , i.e. 
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( 96 ) 



lim (5 (f/(x,t)-F) 


AXy — ^0 


Uf (x + Axj , 0 - [// (jc, ty 
A x j i 


/5(t/(x^)-F)[/ / (jc + Ajc y ,0-5(f/(jc^)-F)[/ i -(jc,0' 


AXy — ^0 \ 


At; 


Now the first term on the R.H.S of (96) involves a two-point correlation between the fine grained PDF at 
the point x and the velocity at the point x + Ax j , which needs a two-point PDF to express the correlation. 
As an approximation (or modeling) for this two-point correlation at the limit A Xj — » 0 , we express it with 
the one-point PDF as follows 

(5 (U(x,t) - V) Uj(x + Axj,tj) = J 5(F' - V) Vif(V';x + Axj,t)dV' (97) 


Then Equation (96) becomes 


y dxj / Ax,--> qJ A Xj 


— lim 

A Xj — ^0 


r Vifu (Vjx + Axj , Q - VJu (V; x, t) ' 
Ax i 


= V t 


dfu(V;x,t) 

dXj 


(98) 


We note that Equation (98) itself is not exact correct due to the assumption made in Equation (97). On the 
other hand, the following integral equation is exactly correct: 



(99) 


This illustrates that the model (98) is just a sufficient condition for Equation (99) to be satisfied. 
Equation (98) also directly suggests 




J_Vu 
fu dxj 


(100) 


Using the same procedure to carry out the mean on the right hand side of (94), we may obtain 



1 dP d 2 U i | 

hi/ — 

p dx j dxjdxj 


a( P(r)df u iv r,a 2 fP 

dV. ^ p dx i dxjdxj y 


(101) 


Therefore, we obtain a closed velocity PDF equation: 

dfu iV dfu _ d ( P(V) dfu + 
dt j dxj dV[ p dxj dxjdxj ^ 


( 102 ) 


Note that the conditional means appeared in the traditional PDF equations (Pope (Ref. 1)) are related to 
the terms on the right hand side of Equation (102). For example, 
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d 2 Uj 

dxjdxj 


-fu 


d 2 Uj 

dxjdxj 



= v 


Vjd 2 fu 

dxjdxj 


This indicates that the conditional mean is 


Similarly, 


v- 


d 2 U t 

dxjdxj 


V ) = v 


V, tfjv 

fu dx j dx j 


/l dP \_1 d[P{V)fu] 

\P dxj I p fu dxi 


(103) 


(104) 


It is clear that these conditional means are actually the functional forms for V 2 Uj and VP in the sample 
space V . Let us denote them as L(Vj) and Pf (V) , then the velocity PDF equation can be simply written 
as 


dfu j y <¥u 

dt 2 dxj 


d 

dK 


P?(V) 


fu + vL(Vi)fu 


Pf(V) = 


1 d[P(V)fu] 
fu Sxi 


L(Vi) = Vi 


V 2 fu 
fu 


(105) 


Now, let us apply the fine grained PDF to the Poison equation of P(x,t ) , which will lead to a functional 
form of P(V ) : 




f 


, d 2 UjUj 
dxjdxj 


(106) 


The operation of the above mean, under a similar assumption made in Equation (97), gives the following 
equation for P(V) 


!WMd w d lk 

p 1 7 dxjdxj 


(107) 


A.3 Transport Equation for Joint Species- Velocity PDF^®(F,\|/; x,t) 

Similar to Equation (94), the transport equation for the fine grained species PDF 
/'(v;*, 0 s 8 (®(*, 0 -v) is 


dt 


df_ 

1 dx j 


c’Vi 


/T 


8 2 <l>, 

dxjdxj 


f'Sf® i0 2 -) 


(108) 


Analogy to the derivation of velocity PDF equation, we may construct the following species PDF 
equation: 
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, V 2 /c/,<t 

i ~ 5 

Ju,o 


^(¥/) -¥/ 


V 2 /t/,<P 



Appendix B. — Miscellaneous Formulations 

In the derivation of PDF equations, the following relationships have been intensively used. The proof 
of Equation (111) has been described by Pope (Ref. 1 ). The same procedure can be followed to prove 
Equation (112). 


f^r(/i/(F;x,OA(F))rfF = 0 

J dV j 

Where A t (V) is a vector function of sample space variable V, and it has a finite mean. 


f s 

t— [/$(¥; x,t)A i (\ i /)]d^ = 0 

fyj 


Where Aj(\\f) is a vector function of sample space variable \|/ , and it has a finite mean. 
Therefore, the following integrations are all zero: 


dVi 


I 8P 

p dxj 


V )fu(V ; x,t)Vj 


dV = 0 


\~^ Vi [ p <y)fu<y-,*,t )]) dv = 0 

p Or j 

J -^r{Vifu(y\ x,o(vV 2 C/,-|f))</F = 0 

v l L(Vj) fu (V; x,t))dV =Q. 


\|/ f /©(v; AT,0(V 2 Oy|\|/^c/v|/ = 0 


I J V/<S/ (ViV • • ¥« )fu,<s> (V, ¥; x i , 0 ) dvdy = 0 


(111) 


( 112 ) 
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